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ABSTRACT 

V445 Pup (2000) is a unique object identified as a helium nova. Color indexes during the outburst are 
consistent with those of free-free emission. We present a free-free emission dominated light curve model of 
V445 Pup on the basis of the optically thick wind theory. Our light curve fitting shows that (1) the white 
dwarf (WD) mass is very massive (Mwo ?J 1-35 Mq), and (2) a half of the accreted matter remains on the 
WD, both of which suggest that the increasing WD mass. Therefore, V445 Pup is a strong candidate of Type 
la supernova progenitor The estimated distance to V445 Pup is now consistent with the recent observational 
suggestions, 3.5 < < 6.5 kpc. A helium star companion is consistent with the brightness of my = 14.5 mag 
just before the outburst, if it is a little bit evolved hot (log T {K) > 4.5) star with the mass of Mee ^ 0.8 Mq. We 
then emphasize importance of observations in the near future quiescent phase after the thick circumstellar dust 
dissipates away, especially its color and magnitude to specify the nature of the companion star. We have also 
calculated helium ignition masses for helium shell flashes against various helium accretion rates and discussed 
the recurrence period of helium novae. 

Subject headings: binaries: close — novae, cataclysmic variables — stars: individual (V445 Pup) — stars: 
mass loss — white dwarfs 



1. INTRODUCTION 

The outburst of V445 Pup was discovered on 30 Decem- 
ber 2000 by Kanatsu (Kato & Kanatsu 2000). The outburst 
shows unique properties such as absence of hydrogen, un- 
usually strong carbon emission lines as well as strong emis- 
sion l i nes of Na, Fe, Ti, Cr, Si, M g etc. jAshok & Banerjeel 
120031: lliiima & NakanishilllOO T'). The spectral features re- 
semble those of classical slow novae exce pt absence of hydro- 
gen a nd strong emission lines of carbon (llijima & Nakanishil 
l2007h . The development of the light curve is very slow (3.3 
mag decline in 7.7 months) with a small outburst amplitude 
of 6 mag. From these features, this object has been suggested 
to be the first example of he lium novae ( Ashok & Banerjee 
HoollKatoAHachisu "2003').^^__^ 

In our earlier work (Kato & Hachisull2003l) we presented 
a theoretical light curve model with the assumption of black- 
body emission from the photosphere, which resulted in the 
best fit model of a very massive white dwarf (WD) (Mwd ^ 
1.3 Mjr,) an d a relatively short distance of c/ < 1 kpc. How- 
ever, lliiima & Nakanishi (2007) recently showed that there 
are absorption lines of Na I Dl/2 at the velocities of 16.0 and 
73.5 km s"', suggesting that V445 Pup is located in or be- 
yond the Orion arm and its distance is as large as 3.5-6.5 
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kpc. IWoudt & Steeghsl (l2008h also suggested a large distance 
of 4.9 kpc. Moreover, color indexes during the outburst are 
consistent with those of free-free emission. With these new 
observational aspects we have revised the light curve model 
of V445 Pup. 

Helium novae were theoretically predicted by Kato et aiTI 
(11989) as a nova outburst caused by a helium shell flash on a 
white dwarf. Such helium novae have long been a theoretical 
object u ntil the first helium nova V445 Pup was discovered 
in 2000. iKato et alTI (1989) assumed two types of helium ac- 
cretion: one is helium accretion from a helium star compan- 
ion, and the other is hydrogen-rich matter accretion with a 
high accretion rate and, beneath the steady hydrogen burning 
shell, ash helium accumulates on the WD. The latter case is 
divided into two kinds of systems depending on whether hy- 
drogen shell burning is steady (stable) or not. To summarize, 
helium accretion onto a WD occurs in three types of systems: 
(1) helium accretion from a helium star companion; (2) steady 
hydrogen accretion with a rate high enough to keep steady hy- 
drogen shell burning s uch as supersoft X-ray sources, e.g., RX 
J05 13 .9-6951 (e.g., Hachisu & Kato 2003; McGowan eLaLl 
l2006h and CAL 83 (e.g., Schmidtke & Cowle y 200 6). (3) hy- 
drogen accretion with a relatively low rate resulting in a re- 
current nova such as RS Oph (e.g ., lHachisu et al. Il2007h and 
U Sco (e.g.. iHachisu et al.1l2000l) . In the present paper, we 
regard that V445 Pup is Case (1), because no hydrogen lines 
were detected (lijima & Najkanishi 200_7). 

In helium novae, mass loss owing to optically thick wind 
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Fig. 1. — V and /c light curves of V445 Pup (lai'ge open circles). The 
/c data are shifted up by one magnitude in order to separate it from the V 
light curve. Open triangles denote V magnitudes (lAUC No. 7552, 7553, 
7557, 7559, 7569, 7620). Small open circles and downward aiTows are taken 
from AAVSO to show very early and later phases of the outburst. Large 
downward arrows show an up per limit ob servation of Ic > 16.2 (this work), 
V >20 and / > 19.5 (Henden e t al. I200TI) . Straight lines indicate the average 
decline rates of 1.9 mag/ 130 days = 0.0146 mag day"' and 3.3 mag/230 days 
= 0.0143 mag day"' for /c and V, respectively. 

is relatively weak compared to hydrogen novae, and a large 
part of the accreted heliu m burns into carbon and oxygen and 
accumulates on the WD (iKato & Hachisu Ill999l [2004 ). Our 
previous model indicated that V445 Pup contains a massive 
WD (MwD ^1-3 Mq) and the WD mass is increasing through 
helium shell flashes. Therefore, V445 Pup is a strong can- 
didate of Type la supernova progenitors. Although there are 
two known evolutional paths of the sing le degenerate scenario 
toward Type la supernova (e.g., Hac hisu. Kato & Nomotol 
[19991) . helium accreting WDs such as V445 Pup do not be- 
long to either of them. Therefore, it may indicate the third 
path to Type la supernovae. If various binary parameters of 
V445 Pup are determined, they provide us important clues to 
binary evolutions to Type la supernovae. 

In the next section we introduce our multi-band photomet- 
ric observations, which indicates that free-free emission dom- 
inates in optical and near infrared. In §3 we present our free- 
free emission dominated light curve model and its application 
to V445 Pup. Then we discuss brightness in quiescence in §4. 
Discussion and conclusions follow in §5 and §6, respectively. 

2. OBSERVATION 

Shortly after the discovery of V445 Pup outburst, one of 
us (S.K.) started multi-band photometry with a 25.4 cm tele- 
scope [focal length = 1600 mm, CCD camera = Apogee AP-7 
(SITE SIA502AB of 512 x 5 12 pixels)]. V andR^ magnitudes 
are obtained from January 4 to May 6, 2001, and 1^ from Jan- 
uary 4, 2001 to January 15, 2007, with the comparison star, 
TYC 6543-2917-1 (V = 8.74, B-V = 0.304). All of our data 
can be taken from the data archive of Variable Star Observers 
League in Japan (VSOLJ)'. 

Figure [T] shows our V and magnitudes as well as other 
observations taken from literature. These two light curves 
show very slow evolution (~ 0.014 mag day"') followed 
by an oscillatory behavior in the V magnitude before it 
quickly darkened by dust blacko ut on about JD 2452 1 00, i.e., 
7.5 months after the discovery dHenden et al. Il200lt IWoudt I 

' http://vsolj.cetus-net.org/ 
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Fig. 2. — Color indexes of V -/c and I' -Sc as well as /c and y. Open 
circles: this work. Open triangles: Gilmore (lAUC 7559, 7569) and Gilmore 
& Kilmartin (lAUC 7620). The horizontal dashed Unes indicate the mean 
values ofV-Rc = 0.5 and V-Ic = 1.12. 



I2002t lAshok & Banerjeell2003h . Here, we assume the de- 
cline rate of the light curves as shown in Figure[T]ignoring the 
later osciflatory phase (JD 2452040 - 2452100) because we 
assume steady-state in the nova envelope as mentioned below 
and, as a result, our theoretical light curve cannot treat un- 
steady oscillations. 

Figure [2] shows evolution of color indexes V -Ic and V-Rc 
as well as V and 1^ themselves. We can see that both ofV-Ic 
and V -Rc are almost constant with time, i.e., « 1.12 and w 
0.5, respectively, during our observation. This means that the 
three light curves of V, /c, and R^ are almost parallel and each 
light curve shape is independent of wavelength. These are the 
characteristic properties of free-free emission dominated light 
curves as explained below. 

The flux of optically thin free-free emission is inversely 
proportional to the square of wavelength, i.e., F\ oc A~^. This 
spectrum shape is unchanged with time although the electron 
temperature rises during nova outburst because the emission 
coefficient dep ends only very slightly on the electron tem- 
perature (e.g.. Brown and Mathews 1970). Therefore, the 
color indexes of free-free emission dominated light curves 
are unchanged with time. On the other hand, the color in- 
dexes change with time if blackbody emission dominates be- 
cause the photospheric temperature rises with time in nova 
outbursts. Figures [1] and |2| indicate that the emission from 
V445 Pup is dominated by free-free emission rather than by 
blackbody emission. 

The color index of free-free emission is calculated from 
XFx oc A~' . When the reddening is known, the reddened color 
is obtained as 



mv-mx = {Mv-Mx\) + cxE(B-V), 



(1) 



where (Mv-Ma)o is the intrinsic color and cx is the red- 
dening coeffic i ent (th ese values are tabulated in Table 4 in 
iHachisuetalTI (l2007h for five colors of V-R, V-I, V-J, 
V-H, and V-K). Figure [3] shows the color indexes relative 
to V for two cases of reddening, i.e., E(B-V) = and 0.5. 

For compariso n, we have added fo ur relative colors 
of VI 500 Cyg (iHachisu&Katol [20061) because its con- 
tinuum flux is known to be that of free-free emission 
dGaUagher & Nevlll976l) except for the first few days after 
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Fig. 3. — Color indexes against wavelength. Our observational color 
indexes of V445 Pup, V -Rc and V —Ic, are plotted for individual data (dots) 
and their mean values of 0.5 and 1.12 (open squares). The color indexes of 
VI 5 00 Cyg are add ed for y magnitude instead of V (filled triangles) (taken 
from lHachisu et alTl ilOO?). Arrows indicate the effective wavelength of six 
bands of V, R, I, J, H, and K. Solid curves denote color indexes of free-free 
emission (or E{B-V) = and 0.5. 

the optical peak. Here we use Stromgren y magnitude in- 
stead of V, because the intermediate-wide y band is almost 
emission-line free. Three {y-I, y-H, and y-K) of four 
color indexes are consistent with those of color indexes with 
£'(B-y) = 0.5 (Ferland 1977; Hachisu & Kato 2006). The 
J band is stron gly contaminated w ith strong emissio n lines 
of He l (e.g.,Black & Gallagher 1 1 976; .Kolotilov & LibermanI 
fT976tlShenavrin & Morozlll976h . This is the reason why the 
V - 7 color index devi ates from that of free-free emission (see 
iHachisu. Kato. & Ca ssatella 2007, for more details). 

Our observed color indexes are shown in Figure |3] The 
open squares are centered at the mean values ofV-Rc = 0.5 
and y-/c = 1.12, whereas dots represent individual obser- 
vations. In V445 Pup, no strong emission lines dominate 
the c o ntinuum in V, R, an d / bands ( Kamath & Anupama] 
I2002t lljjima & Nakanishi 2007). Therefore we conclude 
that the color indexes of V445 Pup are consistent with those 
of free-free emission with the reddening of E{B-V) = 0.51 
(llijima & Nakanishi IIIOOTh . 

3. MODELING OF NOVA LIGHT CURVES 

3.1. Basic Model 

After a thermonuclear runaway sets in on the surface of 
a WD, the envelope expands to a giant size and the opti- 
cal luminosity reaches its maximum. Optically thick winds 
occur and the envelope reaches a steady state. Using the 
same method and numerical techniques as in^ Kato & Hachisul 
([1994, 2003), we have followed evolution of a nova by con- 
necting steady state solutions along the sequence of decreas- 
ing envelope mass. We have solved the equations of mo- 
tion, continuity, radiative diffusion, and conservation of en- 
ergy, from the bottom of the helium envelope through the pho- 
tosphere. The winds are accelerated dee p inside the photo- 
sphere. Updated OPAL opacities are used (llglesias & Rogersi 
11996) . As one of the boundary conditions for our numeral 
code, we assume that photons are emitted at the photo- 



TABLE 1 
Model Parameters" 



MwD 


logRwD 


Y 


Xc+o 


(Mq) 


(«©) 






1.2 


-2.193 


0.68 


0.3 


1.3 


-2.33 


0.48 


0.5 


1.33 


-2.417 


0.38 


0.6 


1.35 


-2.468 


0.38 


0.6 


1.37 


-2.535 


0.58 


0.4 


1.37 


-2.535 


0.38 


0.6 


1.37 


-2.535 


0.18 


0.8 


1.377 


-2.56 


0.38 


0.6 


" The heavy element content is assumed 



to be Z = 0.02 that includes carbon and 
oxygen by the solar ratio. 

sphere as a blackbody w ith a photospheric temperature, Tph. 
iKato & Hachisu I (l2003l) calculated the visual magnitude My 
on the basis of blackbody emission and constructed the light 
curves. In the present work, we calculate free-free emission 
dominated light curves. The envelope structure, wind mass 
loss rate, photospheric temperature, and photospheric radius 
of the WD envelope are essentially the same as those in the 
previous model. 

The flux of free-free emission from the optically thin re- 
gion outside the photosphere dominates in the continuum flux 
of optical and near infrared wavelength region and is approx- 
imated by 

r r°° m2 lull 

cx / NeNidy « / ^^^s^r^dr oc (2) 

J Jr., ^wind'^ V^Ph 

where /v is the flux at the frequency i/, Ne and A^, are 
the number densities of electrons and ions, respectively, V 
is the volume, R^h is the photospheric radius, M„ind is the 
wind mass loss rate, and Vph is the photo spheric velocity 
(IHachisu & Kato Il2006t IHachisu et al. ||2007|) . 

The proportionality constant in equation (|2]i cannot be 
determined a priori because we do not calculate radia- 
tive transfer outside the photosphere: These proportion- 
ahty constant are determined by fitting with observa- 
tional data ( Hachisu & Kato II2006I; IHachisu and Kato II2007I; 
IHachisu et al. Ii2007() . 

3.2. Free-free Light Curve and WD Mass 

We have calculated free-free emission dominated light 
curves of V445 Pup for WD masses of 1.2, 1.3, 1.33, 1.35, 
1.37, and 1.377 Mp, . The last one is the upper limit of the 
mass-accreting WDs (iNomoto et al. 11984) . The adopted WD 
radius and chemical composition are listed in Table [T] Here 
we assume that chemical composition is constant throughout 
the envelope. 

Figure m shows the calculated free-free light curves. More 
massive WDs show faster decline. This is mainly because the 
more massive WDs have the less massive envelope and the 
smaller mass envelope is quickly taken off by the wind. The 
figure also shows two additional models of Mwd = 1-37 Mq 
with different compositions of Xq+o = 0.4 and 0.8. These 
models show almost similar decline rates to that of Mwd = 
1.37 M© with Xc+o = 0.6. 

The starting point of our model light curve depends on the 
initial envelope mass. When an initial envelope mass is given, 
our nova model is located somewhere on the light curve. For 
a more massive ignition mass, it starts from an upper point 
of the light curve. Then the nova moves rightward with the 
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Fig. 4. — Free-free emission dominated liglit curves for vaiious WD 
masses. The cliemical composition is assumed to be 7 = 0.38, Xc+o = 0.6, 
and Z = 0.02 throughout the envelope. The WD mass is attached to each 
curve. The dashed and dotted curves denote Mwd = 1-37 Mq model but with 
different chemical composition of Xq+q = 0.4 and 0.8, respectively^^ Small 
open circles denote the starting points of light curve fitting in Fig. \5\ upper 
and lower points of Mwd = 1 .37 Mq (solid line) correspond to Models 4 and 
3 and those of 1 .377 Mq con'espond to Models 8 and 7, respectively. A short 
solid line indicates the decline rate of observed V data in Fig. [T](3.3 mag / 
230 days). 

decreasing envelope mass due to wind mass loss and nuclear 
burning. 

The development of helium shell flashes is much slower 
than that of hydrogen shell flashes mainly because of 
much larger envelope masses of helium shell flashes 
dKato & Hachisu 111 9941) . The short straight solid line in Fig- 
ure |4]represents the decline rate of V light curve in Figure [T] 
The length of the line indicates the observed period which is 
relatively short because of dust blackout 230 days after the 
optical maximum. We cannot compare the entire evolution 
period of the light curve with the observational data so that 
fitting leaves some ambiguity in choosing the best fit model. 
Even though, we may conclude that WDs with masses of 
^WD ^ 1-33 Mq are very unlikely because their light curves 
are too slow to be comparable with the observation. 

Figure|5]represents the light curve fitting with observational 
data. Here, the model light curves of 4 are identical to those 
of V but are lifted up by 1.12 mag, which is the mean value of 
V-Ic obtained in Figure|2] It should be noted that, in this fig- 
ure, these light curves are further lifted up by 1 .0 mag because 
of/,.-l. 

These light curves are a part of the model light curves in 
Figure |4] In the early phase of the outburst the light curve 
declines almost linearly so that fitting is not unique. We can 
fit any part of our model light curve if its decline rate is the 
same as « 0.014-0.015 mag day"', for example, either top or 
middle part of the same light curve of 1 .377 Mq WD in Figure 
|4] In such a case we show two possible extreme cases, i.e., 
the earliest starting point and the latest one by open circles, as 
shown in Figured These model parameters are summarized 
in TablelH where we distinguish the model by model number 

Thus we have selected several light curves which show a 
reasonable agreement with the observation. However, we can- 
not choose the best fit model among them because the ob- 
served period is too short to discriminate a particular light 
curve from others because the difference among them appears 
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Fig. 5. — Our model light curves are fitted with the observation. Dotted 
line: Model 1. Thin dashed line: Model 2. Thin solid line: Model 3. Thick 
solid line: Model 4. Dash-three dotted line: Model 5. Thin dash-dotted line: 
Model 6. Dashed line: Model 7. Dash-dotted line: Model 8. Our model 
/c light curves are identical with those of V but are lifted up by 1.12 mag 
(and further shifted up by 1 mag for /c — 1). The ordinate on the right vertical 
axis represents the absolute magnitude of Model 4. For the other models, the 
model light curves are shifted down by 0.4 mag (Model 1), up by 0.2 mag 
(Model 2) up by 0.3 mag (Models 3 and 7), down by 0.2 mag (Model 5), 
down by 0.5 mag (Model 6), up by 0.1 mag (Model 8). The observational 
data are same as those in Fig. [T] 

only in a late stage. 

For relatively less massive WDs of Mwd S 1-33 Mq, we 
cannot have reasonable fits with the observation for any part 
of the model light curves. Among relatively more massive 
WDs of Mwd > 1-35 M© in Table El Model 1 shows sHghtly 
slower decline. Therefore, the 1.35 Mq WD may be the low- 
est end for the WD mass. Thus we may conclude that V445 
Pup contains a very massive WD of Mwd ^ 1-35 Mq. 

3.3. Mass Accumulation Efficiency 

During helium nova outbursts, a part of the helium enve- 
lope is blown off in winds, while the rest accumulates on the 
WD. We define the mass accumulation efficiency, yyHe, as the 
ratio of the envelope mass that remains on the WD after the 
helium n ova outburst to the heliu m envelope mass at ignition, 
AMHe.ig (lKato& Hachisu 112004 . 

The mass accumulation efficiency is estimated as follows. 
We have calculated the mass lost by winds during the out- 
burst, AMej. Note that AMej is the calculated total ejecta 
mass which is ejected during the wind phase, and not the mass 
ejected during the observing period which may be shorter than 
the wind phase. 

The ignition mass is approximated by the envelope mass 
at the optical peak, i.e., AMndg ~ AMncpeak- Since the ex- 
act time of the optical peak is unknown, we assume that the 
optical peak is reached on JD 2451872, i.e., the earliest pre- 
discovery observation in the brightest stage reported to lAU 
Circular No. 7553, 28 days earlier than the discovery. The 
envelope mass at the optical peak is estimated from our wind 
solution and is listed as AMHe,ig in Table |2] The resultant 
accumulation efficiency, 

AMHe.peak-AMej 

^He = ITTf , (3) 

^™ He, peak 

is also listed in Table |2] The efficiencies are as high as 
~ 50%. The WD of V445 Pup is already very massive 
(Mwd ^ 1 -35 Mq) and its mass has increased through helium 
nova outbursts. Therefore, V445 Pup is a strong candidate of 
Type la supernova progenitors. 
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TABLE 2 
Parameters of Fitted Model 



model 


MwD 


Y 


Xc+o 






'?Hc''' 


Distance 




' rec 


No. 


(M0) 






(Mq) 


(Mq) 




(kpc) 


(A/oyr-l) 


(yr) 


1 


1.35 


0.38 


0.6 


3.3E-4 


1.8E-4 


0.45 


6.6 


6.E-8 


5000 


2 


1.37 


0.58 


0.4 


1.6E-4 


9.2E-5 


0.42 


5.0 


7.E-8 


2000 


3 


1.37 


0.38 


0.6 


1.9E-4 


8.7E-5 


0.53 


4.8 


6.E-8 


3000 


4 


1.37 


0.38 


0.6 


2.1E-4 


l.lE-4 


0.49 


5.5 


5.E-8 


4000 


5 


1.37 


0.18 


0.8 


3.5E-4 


1.2E-4 


0.64 


6.0 


4.E-8 


9000 


6 


1.37 


0.18 


0.8 


3.8E-4 


1.5E-4 


0.61 


6.9 


4.E-8 


10000 


7 


1.377 


0.38 


0.6 


1.8E-4 


7.4E-5 


0.60 


4.8 


5.E-8 


4000 


8 


1.377 


0.38 


0.6 


2.2E-4 


9.6E-5 


0.56 


5.2 


4.E-8 


5000 



^ Envelope mass at ignition is assumed to be equal to the mass at the optical peak, which we assume as JD 245 1 872 
(28 days before the discoveiy).^ Total mass ejected by the wind'^ Mass accumulation efficiency'^ estimated from 



Fig. El 

4. QUIESCENT PHASE 

Before the 2000 outburst, there was a 14.5 mag star at the 
position of V445 Pup (taken from the archive of VSNET^), 
but no bright star has been observed since the dust blackout. 
We may regard that 14.5 mag is the preoutburst magnitude of 
the binary. There are two possible explanations for this qui- 
escent phase luminosity, one is the accretion disk luminosity 
and the other is the luminosity of bright companion star. In the 
following subsections, we discuss how these possible sources 
contribute to the quiescent luminosity. 

4.1. Accretion Disk 

In some nova systems, an accretion disk mainly contributes 
to the brightness in their quiescent phase. If the preoutburst 
luminosity of V445 Pup comes from an accretion disk, its ab- 
solute magnitude is approximated by 



Mv(obs) = -9.48-- log , 

3 ^ V ^0 yr" 

-2log(2cos/), 



(4) 



where Mwd is the WD mass, Ma^c the mass a ccretion rate, an d 
i the inclination angle (equation (A6) in Webbink et al.ll98"7h. 
Assuming that Mwd = 1 -37 Mq and / = 80° (IWoudt & Steeghsl 
120081) . we have My = 1.4, 2.3, and 3.1 for the accretion rates 
of 1 X 10"'', 3 X 10"^, and 1x10"^ Mq yr"', respectively. 
The apparent magnitude of the disk is calculated from 

m^-Mv = Ay + 5 logDio, (5) 

where Dio is the distance divided by 10 pc. With the absorp- 
tion of Ay = 1 .6 and the distance of 3 kpc, the apparent magni- 
tudes estimated tobemv= 15.4, 16.3 and 17.1, for the above 
accretion rates, respectively. For the distance of 6.5 kpc, we 
obtain brightnesses of my = 17.1, 18.0, and 18.8, respectively. 
All of these values are much fainter than 14.5 mag. Therefore, 
it is very unlikely that an accretion disk mainly contributes to 
the preoutburst luminosity. 

4.2. Helium Star Companion 

Another possible source of the preoutburst brightness is a 
helium star companion. Figure |6] shows evolutional tracks of 
helium stars with masses between 0.6 and 3.0 Mq from the 
helium main-sequence to the red giant stage in the H-R dia- 
gram. We use OPAL opacities and an initial chemical com- 
position ofX = 0.0, Y = 0.98 and Z = 0.02. T he numerica l 
method and input physics are the same as those in lSaiol(ll995h . 

^ http://vsnet.kusastro.kyoto-u.ac.jp/vsnet/ 
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Fig. 6. — Evolutional tracks of helium stars with masses of Mhc = 0.6, 
0.7, 0.8, 0.9, 1.0 (thick line), 1.2, 1.4, 1.6, 1.8, 2.0 (thick line), 2.5, and 
3.0 Mq in the H-R diagram. The 0.6 and 0.7 Mq stars evolve blueward and 
do not become red giants. We stopped calculation when carbon ignites at the 
center for 2.5 and 3.0 Mq stars. Open circles and squares denote stars with a 
14.5 mag brightness at the distance of 6.5 kpc and 4.9 kpc, respectively, for 
Av = 1.6. 

As shown in Figure|6l low mass helium stars do not evolve 
to a red giant. In our new calculation, the 0.8 Mq helium star 
evolves to a red giant, b ut the 0.6 and 0.7 Mq helium stars 
evolve toward blueward. Paczynski I (1197 Ih showed that stars 
of Mile 10 Mq evolve to a red giant while 0.5, 0.7, and 
0.85 Mq do no t. Our calcu lations are essentially the same 
as those of Pacz vrisldl ( 1197 Ih and the difference is attributed 
mainly to the difference between the adopted opacities. 

Figure|6]also shows locations of stars whose apparent mag- 
nitudes are my = 14.5 for Ay = 1.6. Here, the distance is as- 
sumed to be 4.9 (open squares) or 6.5 kpc (open circles). For 
example, a star of (temperature, luminosity) = (log Tph (K), 
logLph/LQ) = (4.2, 2.81), (4.4, 3.25), and (4.6, 3.75) could 
be observed as a 14.5 mag star for 6.5 kpc and (logTph (K), 
logLph/LQ) = (4.2, 2.57), (4.4, 3.0), and (4.6, 3.51) for 4.9 
kpc. Therefore, the observed 14.5 mag is consistent with 
the luminosities and temperatures of slightly evolved helium 
stars of Mhb > 0.8 Mq if the companion is a blue star of 
log7ph>4.5. 

On the other hand, if the companion is redder than log Tph < 
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Fig. 7. — The helium ignition mass, AMncig, of helium-accreting WDs is 
plotted against the helium accretion rate, Mhc ■ The WD mass is attached to 
each curve. Straight dashed lines indicates the recurrenceneriod. The open 
square indicates the ignition mass of each model in Table|2] 

4.4, there is no corresponding evolution track in the low lumi- 
nosity region of the H-R diagram as shown in Figure |6] In 
such a case the preoutburst luminosity cannot be attributed to 
a helium star companion. Kato (2001)-' suggested that preout- 
burst color of V445 Pup was much bluer than that of symbi- 
otic stars (which have a red giant companion) but rather close 
to that of cataclysmic variables (main-sequence or little bit 
evolved companion). This argument is consistent with our es- 
timate of relatively high surface temperature log Tph > 4.5 (see 
Fig.©. 

5. DISCUSSIONS 
5.1. Distance 

As explained in §§3.1 and 3.2, we do not solve energy trans- 
fer outside the photosphere and, thus, we cannot determine the 
proportionality constant in equation (|2]l. Therefore, we cannot 
determine the distance to V445 Pup directly from the compar- 
ison of theoretical absolute magnitude with the observational 
apparent magnitude. Instead, we can approximately estimate 
the distance by assuming that the free-free flux is larger than 
the blackbody flux during the observing period. This assump- 
tion gives a lower limit to the distance. We expect that the dis- 
tance estimated thus is close to the real value. The 8th column 
of Table |2] lists the distance estimated for each model. These 
distances are consistent with the o bservational estimates of 
3.5 < J<6.5kp c (liiim a & Nakanishi 2007) and ^ 4.9 kpc 
jWoudt & Steeghs II2008I) . 

5.2. Helium Ignition Mass and Recurrence Period 

We have calculated evolution of C-i-O white dwarfs accret- 
ing helium at various rates until the ignition of a helium shell 
flash. Chemical compositions assumed are Xc = 0.48 and 
Xo = 0.50 for the WD core, and X = 0.0 and Y = 0.98 for the 
accreted envelope. The initial model adopted for a given core 
mass and accretion rate is a steady-state model in which the 
heating due to the accretion is balanced with the radiative en- 
ergy flow. Adopting such an initial model is justified because 

' http://www.kusastro.kyoto-u.ac.jp/vsnet/Mail/alert5000/msg00493.html 



the WD has accreted matter from the companion for a long 
time and experienced many shell flashes. 

As the helium accretion proceeds, the temperature at the 
bottom of the envelope gradually increases. When the temper- 
ature becomes high enough the triple-alpha reaction causes a 
shell flash. The mass of the helium envelope at the ignition 
depends on the WD mass and the accretion rate as shown in 
Figure |2l The envelope mass required to ignite a shell flash 
tends to be smaller for a higher accretion rate and for a larger 
WD mass. Helium flashes are weaker for higher accretion 
rates and lower WD masses. In particular they are very weak 
for accretion rates higher than ~ 10~^Mq yr"' . 

The envelope mass at the ignition depends also on the core 
temperature. If the core temperature of the initial WD is lower 
than that of the steady state model employed for our calcula- 
tions, the envelope mass at the ignition would be larger. 

The recurrence period, AMne.ig/MHe, is also plotted in Fig- 
ure [T] by dashed lines. The recurrence period corresponding 
to the ignition mass of our fitted models are listed in Table|2l 

5.3. Mass Transfer from Helium Star Companion 

We have estimated the mean accretion rate Mhc of the WD 
as in TablelH The companion star feeds its mass to the WD via 
the Roche lobe overflow or by winds depending on whether it 
fills the Roche lobe or not. Here we examine if this accretion 
rate is comparable to the mass transfer rates of the Roche lobe 
overflow from a helium star companion. We followed stellar 
evolutions of helium stars from the main sequence assuming 
that the star always fills its effective Roche lobe of which ra- 
dius is simply assumed to be 1.5 Rq. The resultant mass loss 
rates are shown in Figure[8]for stars with initial masses of 1 .2, 
1.0 and 0.8 Mq. These results are hardly affected even if we 
change the Roche lobe radius. 

The mass loss rate decreases with the companion mass al- 
most independently of the initial value. Except for the short 
initial and final stages, these rates are as large as Mhb ~ 
10"^ Mq yr"' and much larger than the mass accretion rate 
estimated from our fitted models in Table ID Since almost all 
of the mass lost by the companion accretes onto the WD in 
the case of Roche lobe overflow, the accretion may result in 
very weak shell flashes for such high rates as 1O~''M0 yr"'. 
A strong shell flash like V445 Pup is realized only when the 
final stages of mass-transfer in which the mass transfer rate 
quickly drops with time. This rare case may correspond to a 
final stage of binary evolution after the WD had grown up in 
mass to near the Chandrasekhar limit from a less massive one 
and the companion had lost a large amount of mass via Roche 
lobe overflow. Such information will be useful for modeling a 
new path of binary evolution scenario that will lead to a Type 
la supernova. 

5.4. Comparison with Other Observation 

The mass of the dust shell can be estimated from in- 
frared continuum flux with the assumption that the emis- 
sion o riginated from warm dust. With infrared 10 jim spec- 
trum. Lynch et al. I (L2004b) estimated the dust mass to be 
2 X 10"^Mo for a distance of 3 kpc. This value would be 
increased if we adopt a larger distance instead of 3 kpc, 
but still consistent with our estimated ejecta mass AMgj ^ 
(0.7-1.8) X \0-^Mq in Table |2l because the dust mass IS a 
sr nall part of the ejecte d mass. 

iLvnch et al. I (l2004ah reported the absence of He II or coro- 
nal lines on their near IR spectra and suggested that the ion- 
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Fig. 8. — Mass loss rate of a helium stai' that fills its Roche lobe of the 
radius 1.5 Rq. The initial stellar mass is 1.2, 1.0 (dotted line), and 0.8 Mq. 
Time goes on from left to right. See the text for more detail. 

izing source was not hot enough on 200 4 January 14 and 16 
(1 146 days after the optical maximum). iLynch et al. I (l2005h 
also reported, from their near IR observation on 2005 Novem- 
ber 16 that the object had faded and the thermal dust emission 
had virtually disappeared. They suggested that the dust had 
cooled significantly until that date (1818 days). This sugges- 
tion may constrain the WD mass of V445 Pup because Fig- 
ure |4] shows the outburst lasts more than 1800 days for less 
massive WDs (< 1.33 Mq). If the WD had cooled down un- 
til the above date we may exclude less massive WD models 
(< 1.33 Mq) because these WDs evolve slowly and their hot 
surfaces emit high energy photons at least until 1800 days af- 
ter the optical peak. 

We have shown that the WD mass of V445 Pup is increas- 
ing with time. When the WD continues to grow up to the 
Chandrasekhar limit, central carbon ignition triggers a Type 
la supernova explosion if the WD consists of carbon and oxy- 
gen. We regards the star as a CO WD instead of an O-Ne-Mg 
WD because no indic ation of neon were obs erved in the neb- 
ula phase spectrum (Woudt & Steeghs 2005). Therefore, we 
consider that V445 Pup is a strong candidate of Type la su- 
pernova. 

When a binary consisting of a massive WD and a helium 
star like V445 Pup becomes a Type la supernova, its spec- 
trum may more or less show a sign of helium. The search for 
such helium associated with a Type la supern ova has been re- 
ported for a dozen Type la supernovae (M arion et al. ll200Tt 
iMattila et al 1120051) . but all of them are negative detection. It 
may be difficult to find such a system because this type of 
Type la supernovae are very rare. 

The binary seems to be still deeply obscured by the op- 
tically thick dust shell even several years after the outburst. 



This blackout period is much longer than that of the classical 
novae such as OS And (~ 20 days), V705 Cas 100 days), 
and DQ Her (~ 100 days). As the ejected mass estimated in 
Table |2] is not much different from those of classical novae, 
the difference of the blackout period may be attributed to a 
large amount of dust in the extremely carbon-rich ejecta of a 
helium nova. Moreover, observed low ej ection velocities of 
- 500 km s-i dn] ima & Nakanishi II2007I) is not much larger 
than the escape velocity of the binary with a relatively mas- 
sive total mass [e.g., 1.35 Mo + (l-2) Mq]. Both of these 
effects lengthen the dust blackout period. When the dust ob- 
scuration will be cleared in the future the period of blackout 
provides useful information on the dust shell. 



6. CONCLUSIONS 

Our main results are summarized as follows: 

1 . We have reproduced V and light curves using free-free 
emission dominated light curves calculated on the basis of the 
optically thick wind theory. 

2. From the light curve analysis we have estimated the WD 
mass to be as massive as Mwd ^ 1-35 Mq. 

3. Our light curve models ar e now consistent with a lo nger 
distance of 3.5 <d<6 .5 kpc ( Ijiima & Nakanishi Il2007h and 
4.9 kpc dWoudt & Steeghs .20m . 

4. We have estimated the ejecta mass as the mass lost by 
optically thick winds, i.e., AM^ind ~ 1Q~'^Mq. This amounts 
to about a half of the accreted helium matter so that the accu- 
mulation efficiency reaches ^ 50%. 

5. The white dwarf is already very close to the Chan- 
drasekhar mass, i.e., Mwd ?J 1 -35 Mq, and the WD mass had 
increased after the helium nova. Therefore, V445 Pup is a 
strong candidate of Type la supernova progenitors. 

6. We emphasize importance of observations after the dense 
dust shell will disappear, especially observations of the color 
and magnitude, orbital period, and inclination angle of the 
orbit. These are important to specify companion nature. 
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